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In 2 + 1 flavour lattice QCD the spectrum of the nucleon is presented for both parities using local 
meson-baryon type interpolating fields in addition to the standard three-quark nucleon interpola¬ 
tors. The role of local five-quark operators in extracting the nucleon excited state spectrum via 
correlation matrix techniques is explored on dynamical gauge fields with m n = 293 MeV, leading 
to the observation of a state in the region of the non-interacting S-wave Nn scattering threshold 
in the negative-parity sector. Furthermore, the robustness of the variational technique is exam¬ 
ined by studying the spectrum on a variety of operator bases. Fitting a single-state ansatz to 
the eigenstate-projected correlators provides robust energies for the low-lying spectrum that are 
essentially invariant despite being extracted from qualitatively different bases. 

PACS numbers: 11.15.Ha,12.38.-t,12.38.Gc 


I. INTRODUCTION 

Lattice QCD is currently the only known ab-initio non- 
perturbative approach to study the fundamental quan¬ 
tum field theory governing hadron properties, Quantum 
Chromodynamics (QCD). While the ability to obtain 
ground state masses is well-understood, an accurate ex¬ 
traction of excited states and multi-particle thresholds 
remains a challenge. 

The use of variational techniques mm to study the 
nucleon excited state spectrum has seen remarkable suc¬ 
cess in recent years. The key feature of these techniques 
is to begin with a basis of different operators that couple 
to the quantum numbers of a given state, and then con¬ 
struct different linear combinations of these operators in 
order to isolate the ground and higher excited states in 
that channel. 

The positive-parity nucleon channel has been of signifi¬ 
cant interest to the lattice community SHE]- In particular 
the first positive-parity J p = excitation of the nu¬ 
cleon, known as the Roper resonance TV*(1440), remains 
a puzzle. In constituent quark models the Roper reso¬ 
nance lies above the lowest-lying negative-parity state 0 - 
HTj . the IV* (1535), whereas in Nature it lies 95 MeV 
below the resonant state. This has led to speculation 
about the true nature of this state, with suggestions it 
is a baryon with explicitly excited gluon fields, or that 
it can be understood with meson-baryon dynamics via a 
meson-exchange model [la¬ 
in simple quark models, the Roper is identified with 
an N = 2 radial excitation of the nucleon. Within the 
variational technique, the choice of an appropriate oper¬ 
ator basis is critical to obtaining the complete spectrum 
of low-lying excited states. Recall that we can expand 
any radial function using a basis of Gaussians of differ¬ 
ent widths /(|r|) = Y^i c t e ~ eir • This leads to the use 
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of Gaussian-smeared fermion sources with a variety of 
widths 13], providing an operator basis that is highly 
suited to accessing radial excitations. The CSSM lat¬ 
tice collaboration has used this technique to study the 
nucleon excited state spectrum m ITS) . In particular, 
the CSSM studies were the first to demonstrate that the 
inclusion of very wide quark fields (formed with large 
amounts of Gaussian smearing) is critical to isolating the 
first positive-parity nucleon excited state (8j US]. This 
state was shown to have a quark probability distribution 
consistent with an N = 2 radial excitation in Ref. HZ]. 
This work also examined the quark probability distri¬ 
butions for higher positive-parity nucleon excited states, 
revealing that the combination of Gaussian sources of dif¬ 
ferent widths allows for the formation of the nodal struc¬ 
tures that characterise the different radial excitations. 

The negative-parity nucleon channel with its two low- 
lying resonances, the iV*(1535) and 7V*(1650), has also 
been of significant interest i6] [21 ' 18ff20j . These S'n 
states are in agreement with SU (6) based quark model 
predictions, making an ab-initio study of the low-lying 
negative-parity spectrum a potentially rewarding endeav¬ 
our. Importantly, at near physical quark masses the non¬ 
interacting nN scattering threshold lies below the lowest 
lying negative-parity state, making it a natural place to 
look for the presence of multi-particle energy levels in the 
extracted spectrum. 

Until recently, the majority of the work in these chan¬ 
nels has been performed with three-quark interpolating 
fields, and in the full quantum field theory these interpo¬ 
lators couple to more exotic meson-baryon components 
such as the aforementioned nN via sea-quark loop in¬ 
teractions. However, baryon studies have found that 
the couplings of single hadron type operators to hadron- 
hadron type components, suppressed by the lattice vol¬ 
ume as l/y/V, are sufficiently low so as to make it difficult 
to observe states associated with scattering thresholds 
mm- Moreover, there is a question as to what extent 
the presence of multi-particle states might interfere with 
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the extraction of nearby resonances. 

One solution is to explicitly include hadron-hadron 
type interpolators pa im by combining single-hadron 
operators with the relevant momentum. This creates an 
operator that necessarily has a high overlap with the 
scattering state of interest thereby enabling its extrac¬ 
tion. Instead, in this work we aim to construct meson- 
baryon type interpolators without explicitly projecting 
single-hadron momenta, and investigate the role that the 
resulting operator plays in the calculation of the nucleon 
spectrum. Using these operators we construct a basis 
containing both three- and five-quark operators, and per¬ 
form spectroscopic calculations utilising a variety of dif¬ 
ferent sub-bases. Examining the resulting spectra then 
provides an excellent opportunity to both study the role 
of our multi-particle operators and test the robustness of 
the variational techniques employed. 

Following the outline of standard variational analyses 
in Section [TI| we construct these hadron-hadron type in¬ 
terpolators in the form of five-quark operators in Sec¬ 
tion |III| We then develop a method for smearing ele¬ 
ments of the stochastically estimated loop propagators at 
x, S(x, x) in Section IV These necessarily arise with the 
introduction of our five-quark interpolating fields, due to 
the presence of creation quark fields in our annihilation 
operator and vice versa. Having covered the technology 
required for a spectroscopic calculation, we then outline 


our simulation details in Section [V 
spectra for both parities in Section 


and present nucleon 
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II. CORRELATION MATRIX TECHNIQUES 


Correlation matrix techniques m m are now well- 
established as a method for studying the excited state 
hadron spectrum. The underlying principle is to begin 
with a sufficiently large basis of N operators (so as to 
span the space of the states of interest within the spec¬ 
trum) and construct an N x TV matrix of cross correlation 
functions, 

Gij{p,t) = ^2e~ lp ' x (Cl\xi(x,t)Xj(0,t src )\Cl). (1) 

X 


After selecting p = 0 and projecting to a specific parity 
with the operator 


r± = \ (7o ± I ), (2) 

we can write the correlator as a sum of exponentials, 

= (3) 

a 


respectively. We then search for a linear combination of 
operators 


4> a = Xj Uj and (j) a = Xi v? (4) 

such that (j) and </> couple to a single energy eigenstate. 
That is, we require 


( Q | </>“ |/3) oc <5 a/3 . (5) 

One can then see from Eq. ([3]) that 

Qij{to + dt) uf = e ~ madt Gij(t 0 ) Uj, ( 6 ) 

and hence the required values for and vf for a given 
choice of variational parameters ( to,dt ) can be obtained 
by solving the eigenvalue equations 

[£ -1 (*o)S(*o + di)]..u? = ( 7 ) 

vf [5(t 0 + dt) £ _1 (to)]^. = c a vf, (8) 


where the eigenvalue is c“ = e~ madt . In the ensemble 
average Gij is a symmetric matrix. We work with the 
improved estimator \ (G t j +G ri ) ensuring the eigenvalues 
of Eqs. 0 and Q are equal. As our correlation matrix 
is diagonalised at t 0 and t 0 + dt by the eigenvectors 
and uf we can obtain the eigenstate-projected correlator 
as a function of Euclidean time 

0“(t)=i;?Sy (*)«“> (9) 

which can then be used to extract masses. Moreover, the 
analysis can be performed on a symmetric matrix with 
orthogonal eigenvectors. More details can be found in 
Ref. [ 151 . 

At this point we note that if the operator basis does not 
appropriately span the low-lying spectrum, G a (t) may 
contain a mixture of two or more energy eigenstates. 
There are a number of scenarios in which this might oc¬ 
cur: 

• At early Euclidean times the number of states 
strongly contributing to the correlation matrix may 
be (much) larger than the number of operators in 
the basis. 

• There may be energy eigenstates present that do 
not couple or only couple weakly to the opera¬ 
tors used. In particular, it is well known that lo¬ 
cal three-quark interpolating fields couple poorly to 
multi-hadron scattering states. 

• The nature of the operators selected may be such 
that it is not possible to construct a linear combi¬ 
nation with the appropriate structure to isolate a 
particular state. 


where a enumerates the energy eigenstates of mass m a 
and X° and A“ are the couplings of our creation and 
annihilation operators Xj and Xi a f the source and sink 


It is important to have a strategy to ensure that one can 
accurately obtain eigenstate energies from the correlation 
matrix. The method we use is to analyse the effective 
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energies of different states from the eigenstate-projected 
correlators, 


E a {t ) = - log (10) 

w n & Q a {t + n) K ’ 

which is constant in regions where the correlator is dom¬ 
inated by a single state. Neighbouring time slices in the 
correlation functions are highly correlated in Euclidean 
time, and require a covariance-matrix based y 2 analysis. 
The best unbiased estimate corresponds to a y 2 /dof ~ 1. 
We therefore endeavour to obtain a plateau fit of the ef¬ 
fective mass with the \ 2 /dof close to one. In considering 
an upper limit for the fit, points with errors bars larger 
than the central value are discarded. Fits with y 2 /dof 
> 1.20 are rejected, as these fits have significant contam¬ 
ination from nearby states not yet isolated in the corre¬ 
lation matrix analysis [20]. We do not enforce a lower 
bound on acceptable y 2 /dof as small values typically re¬ 
flect large uncertainties rather than an incorrect result 
associated with a systematic error. Typically, plateaus 
commence three or four time slices after the source, near 
the regime where the generalised eigenvector analysis of 
the correlation matrix is done. Figure[l]illustrates typical 
effective mass fits for positive and negative parity states. 
Further details of this method can found in Ref. [22]. 

As we will demonstrate, a careful covariance-matrix 
based y 2 analysis to fit the single-state ansatz ensures a 
robust extraction of the eigenstate energies. The physics 
underpinning this robustness is elucidated in detail in 
Section IVII 

The CSSM lattice collaboration has used this tech¬ 
nique in the calculation of the nucleon spectra in both the 
positive [5] and negative-parity channels m with stan¬ 
dard three-quark interpolators. While largely successful 
at identifying towers of excited states that would be asso¬ 
ciated with resonances in Nature, it has been shown that 
with three-quark operators alone it is difficult to detect 
states near multi-particle scattering energy levels [20] . 
The concern is that the operator basis doesn’t have suffi¬ 
cient overlap with meson-baryon type components, high¬ 
lighting the need for studies with multi-hadron operators. 


III. MULTI-PARTICLE STATE 
CONTRIBUTIONS 
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FIG. 1. (Colour online). Typical effective mass fits for 
positive-parity (top) and negative-parity (bottom) nucleon 
excitations. The top plot shows a fit to the first positive- 
parity excitation of the 4x4 correlation matrix obtained from 
basis 1 described in Table Q] of Section ED The fitted mass of 
2.11(4) GeV is shown as a green circle for basis 1 in Fig. [4] 
and provides y 2 /dof = 0.17. The bottom plot shows a fit 
to the lowest-lying state in the negative-parity sector. It is 
sourced from the 6x6 correlation matrix obtained from basis 
3 described in Table[I]of Sectionjv] The fitted mass of 1.58(3) 
GeV is shown as a blue square for basis 3 in Fig. [8] and corre¬ 
sponds to y 2 /dof = 0.87. Note, an earlier fit including t = 20 
provides y 2 /dof = 1.22, reflecting the systematic drift in the 
effective mass at early times. 



In order to further elucidate the situation, we consider 
a simple two-component toy model which consists of two 
QCD energy eigenstates, | a) and | b). We then suppose 
that | a) and | b) are given by 

| a) = cos# | l) T sin# | 2) , (11) 

| b) = — sin 6 | l) + cos 0 | 2) , (12) 

where 11) and | 2) denote a single-hadron and meson- 
baryon type component respectively, while # is some ar¬ 
bitrary mixing angle. Now imagine performing a spec¬ 
troscopic calculation with an interpolating field y 3 that 


only has substantial overlap with 11). That is, 

(fl | y 3 I 1 ) CX C and (Q | y 3 | 2) <C C , (13) 

for some constant C. When y 3 acts on the vacuum we 
therefore create a state that is superposition of the true 
energy eigenstates given by 

| l) = cos 6 | a) — sin 9 | b) . (14) 

In the absence of an operator that has substantial over¬ 
lap with | 2), it becomes impossible to separate out the 
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true QCD eigenstates of interest. This naturally leads to 
two points of concern. Firstly, one cannot extract states 
with a significant | 2) component and secondly there is 
possibly contamination of the states that are extracted. 
When performing baryon spectroscopy it therefore be¬ 
comes desirable to include interpolating fields that we 
expect to have substantial overlap with multi-particle 
meson-baryon type states m- While projecting single¬ 
hadron momenta in a multi-hadron operator allows for 
a clean extraction of states associated with scattering 
thresholds, the influence of local five-quark operators 
(without explicit momenta assigned to each hadron) on 
the spectrum is less intuitive. It is the purpose of this 
study to examine the role local five-quark operators play 
in the spectrum, and to thereby test the robustness of 
our variational method. 

Starting with standard N and 7 r interpolators we use 
the Clebsch-Gordan coefficients to project isospin I = 
1/2,1^ = +1/2 and write down the general form of our 
meson-baryon interpolating fields [531 121] , 

XnA x ) = -^e obc 75X 

2\u Ta (x) Ti d 6 (a;)] T 2 d c (x ) \d e (x) 75 u e (x)] 

- [u Ta (x) Ti d b (a;)] r 2 u c (:r) [d e (x) 75 d e (x) ] 

+ [u Ta (x)r 1 d b (x)] r 2 u c (x) [u(x) e ^^{x)] j, 

(15) 

providing us with two five-quark operators, denoted xs 
and X 5 which correspond to (Fi,F 2 ) = (C 7 s,I) and 
(ri,r 2 ) = (C, 75) respectively. The square brackets 
around the diquark contraction denote a Dirac scalar. 
Under a parity transformation 

x — > x = (xq, —x) , (16) 

and the quark fields il>{x) and i/(x) transform as 
ijj(x) -A V tf>(x) V t = 70 ip(x ), 

^j(x) —> V '4){x)V^ = 4>(x)”/o- (17) 



FIG. 2. The Feynman diagrams considered following the in¬ 
troduction of five-quark interpolating fields to standard three- 
quark operators. 

combined in a correlation matrix. This also ensures the 
standard parity projector of Eq. ©> applies to our five- 
quark interpolators. 

The presence of creation quark fields in our annihila¬ 
tion interpolating held and vice versa then leads to the 
requirement of calculating the more computationally in¬ 
tense loop propagators, in order to compute the diagrams 
in Fig. [2] The literature contains different ways of dealing 
with such diagrams such as distillation [25], and various 
schemes such as the Laplacian Heaviside (LapH) smear¬ 
ing method [55] . Here we will stochastically estimate 
inverse matrix elements fully diluting in spin, colour and 
time as outlined below. 


IV. LOOP PROPAGATOR TECHNIQUES 

As observed in the preceding section, spectroscopic cal¬ 
culations that involve the five-quark operators \5 and x's 
necessarily involve the determination of loop propagators 
at x, denoted S(x, x). As S(x, x) requires a source at each 
lattice point, a different recipe to that of conventional 
point-to-all propagators is utilised. For this purpose we 
use stochastic estimation of the matrix inverse [23 [28. 

Given a set of random noise vectors { 77 } with elements 
drawn from Z 4 such that the average over noise vectors 
gives 


Applying a parity transformation to the standard pion 
interpolator Xn{x) = ^(x^si/j^x), and the nucleon inter¬ 
polators of type Xn(x) = [ip T {x)(C^ 5 )ip(x)]4’{x) of Eq. 
(27| we find 


Xn(x) -A -ip(x)'y 5 i/>(x) = -Xtt(x) , 

Xn(x) -a [ip T (x) (C'y 5 )ip(x)]'y 0 ip(x) = 70 Xn(x). (18) 


Thus the pion interpolator transforms negatively under 
parity. To ensure our five-quark baryon interpolator 
formed from the product of pion and nucleon interpo¬ 
lators, transforms in the appropriate manner, the pref¬ 
actor of 75 is included in Eq. (15). That is, both our 
three-quark and five-quark nucleon operators have the 
same parity transformation properties and hence can be 


(jlaot (*£) bb/3 (?/)) — dxydabdaf 3? (1^) 

with colour indices a, b, spin indices a, /3 and space-time 
indices x , y. We define for each noise vector a correspond¬ 
ing solution vector 

X = M^ 1 77 , ( 20 ) 

where in this case M is the fermion matrix. Then the 
stochastic estimate of a propagator matrix element is cal¬ 
culated as 


^ab;ayi,X,y) — ( Xaa(x )dlpiy))- ( 21 ) 

We perform full dilution in time, spin and colour indices 
as a means of variance reduction [5i?] . That is, given a 
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set of full noise vectors {??}, we can define a set of diluted 
noise vectors {rj^ a “ ( 1 } by 

Vaa a ‘ 1 (®) t) = <W <W 8 tt > Vaa (x, t ), (22) 

where the intrinsic quark field indices are specified by 
colour a, spin a, space x and time t respectively and 
the colour-spin-time diluted noise vectors are enumerated 
by the corresponding [o' a! t'\ labels. We can similarly 
enumerate the solution vectors 

x [°'“'t'] =Af—V 0 ' 0 '*'], (23) 

which makes it clear that by diluting we increase the 
number of inversions required by a factor of n C 0 i 0 ur x 
TTgpin x 77-time- The stochastic estimate of the matrix in¬ 
verse with dilution is given by 

S(x,y)^(Y, X la,a,t ' ] (x)v la ' a,t ' ] Hy )>, ( 24 ) 

a' ,ol' ,t' 

where colour and spin indices are taken to be implicit 
for clarity. At this point we remark that, while it is 
computationally infeasible to also fully dilute in the space 
index x, in this extreme limit each diluted noise vector 
would consist of only a single non-zero element, meaning 
that we are exactly calculating the full matrix S(x, y ) and 
the above relation becomes an equality rather than an 
estimate. This makes it clear that using dilution provides 
an improved stochastic estimate to the matrix inverse. 

As shown in Fig. [2} our construction of nucleon correla¬ 
tors with five-quark operators combines standard point- 
to-all propagators S(x, 0) and stochastic estimates of the 
loop propagators S(x,x). In order to access the radial 
excitations of the nucleon, we make use of multiple lev¬ 
els of Gaussian smearing in our quark fields. Hence, to 
construct a correlation matrix we need to calculate prop¬ 
agators with differing levels of source and sink smearing. 

Let (x, y) denote a propagator with m iterations 

of smearing applied at the sink and n iterations applied 
at the source. In the case of point-to-all propagators 
5( m,n )(x, 0 ) the source point is fixed, y = 0 , and starting 
with a point source ip(°\ we apply n iterations of Gaus¬ 
sian smearing pre-inversion to obtain the smeared source 
where 


3 

H i/j(x) = (1 - a) ip(x) + ^ X ^ + a A) 

1 

+U^(x — afi) ip(x — afi)} , (25) 

and a specifies the smearing fraction. Sink smearing 
is applied to the propagator post-inversion to obtain 
S'( m ’")(x,0). 

The application of smearing to construct a stochastic 
estimate for the quark propagator g( m ’ n )(# ? y) is some¬ 
what different. The set of (diluted) noise and solution 
vectors { 77 , x} is first constructed, whereby it follows from 
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FIG. 3. (Colour online). A comparison of correlators calcu¬ 
lated with one stochastically estimated propagator (denoted 
“stochastic”) to those calculated with no stochastic propaga¬ 
tors (denoted “standard”). Results are presented for the pion 
(top) and the ground state nucleon (bottom). 


Eqs. (21) and (251 that an estimate of the smeared prop¬ 
agator is given by 


SM(x,y) = ( X ( - m \x)r 1 ^(y)), (26) 


where x ^ = H m X is the result of m iterations of Gaus¬ 
sian smearing applied to the (diluted) solution vectors, 
and 7 ji"i = H n r] is similarly constructed from the (di¬ 
luted) noise vectors. Note that the smearing is applied 
after (any dilution and) the solution vectors have been 
calculated. The construction of a smeared loop prop¬ 
agator (x, x) is simply an application of the above 

formulae in the case y = x. 

In order to determine how many noise vectors per con¬ 
figuration are sufficient to provide similar statistical er¬ 
rors for our point-to-all and stochastic propagators, cor¬ 
relators are calculated with a stochastic estimate of the 
point-to-all propagator, and compared to those obtained 
using point-to-all propagators calculated in the standard 
way. As each independent quark line in a hadron correla¬ 
tor requires an independent noise source to ensure unbi- 
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ased estimation |26| we insert one stochastic propagator 
into the aforementioned correlators. Furthermore, as a 
test of our smearing technique for stochastic propagators, 
we perform this comparison using a variety of smearing 
levels. Note that, as smearing of both the source and so¬ 
lution vectors is performed post-inversion, the stochastic 
method effectively provides different smearing levels for 
free. 

The comparison is performed on 75 20 3 x 40 lattice 
gauge configurations, with the FLIC fermion action m- 
The lattice spacing is 0.126 fm in both the temporal and 
spatial direction providing a physical lattice volume of 
(2.52 fm) 3 . Four full noise vectors are used per stochas¬ 
tic propagator, which are then colour-spin-time diluted. 
As the source timeslice is fixed in this case, each stochas¬ 
tic propagator requires n co i our x n sp i n inversions per noise 
vector. Recall standard point-to-all propagators require 
^colour x n S pin inversions, although source smearing is 
applied pre-inversion unlike the stochastic case. Three 
different levels of smearing are used, n s = 35,100, 200 
sweeps with a = 0.7. Fig. [3] shows good agreement 
across all smearing levels between those correlators con¬ 
taining a stochastically estimated propagator and those 
that do not, demonstrating that using four noise vectors 
per quark line provides a comparable statistical uncer¬ 
tainty to that of a standard propagator. We note here 
that ultimately we utilise this method to calculate S(x, x) 
not S(x,x source ), meaning we get the added benefit of 
spatial averaging for our loop propagators. 


operators 

Xl = e abc [u aT (C'y 5 )d b ]u c 

X2 = e abc [u aT (C)d b ] l5 u c (27) 

in order to form the seven bases we study that are out¬ 
lined in Table ID 


TABLE I. Table of the various operators used in each basis. 


Basis Number 

Operators Used 

1 

Xl, X2 

2 

XU X2, X5 

3 

Xl, X2, Xs 

4 

XU X2, X5, X5 

5 

Xl, X5, X5 

6 

X 2 , X5, Xs 

7 

xs, Xs 


Throughout this work we employ Gauge-invariant 
Gaussian smearing 1361 at the source and sink to increase 
the basis size via altering the overlap of our operators 
with the states of interest. We choose n s = 35 and 
n s = 200 sweeps of smearing providing bases of sizes 4, 
6 and 8 . Stochastic quark lines are calculated using four 
random Z 4 noise vectors that are fully diluted in colour, 
spin and time. 


V. SIMULATION DETAILS 


VI. RESULTS 


For the baryon spectroscopy results presented herein 
we use the PACS-CS 2+1 flavour dynamical-fermion con¬ 
figurations m made available through the ILDG 32]- 
These configurations use the non-perturbatively O(a)- 
improved Wilson fermion action and the Iwasaki gauge 
action [32!. The lattice size is 32 3 x 64 with a lat¬ 
tice spacing of 0.0907 fm providing a physical volume 
of ss (2.90 fm) 3 . /3 = 1.90, the light quark mass is set 
by the hopping parameter n u d = 0.13770 which gives a 
pion mass of m n = 293 MeV, while the strange quark 
mass is set by k s = 0.13640. Fixed boundary conditions 
are employed in the time direction removing backward 
propagating states [311 S 3 ] , and the source is inserted at 
t src — n t /4 = 16, well away from the boundary. System¬ 
atic effects associated with this boundary condition are 
negligible for t > 16 slices from the boundary. The main 
results of our variational analysis is performed at to = 17 
and dt = 3, providing a good balance between systematic 
and statistical uncertainties. Uncertainties are obtained 
via single elimination jackknife while a full covariance 
matrix analysis provides the x 2 /dof which is utilised to 
select fit regions for the eigenstate-projected correlators. 


In addition to the five quark operators \s and X 5 pre¬ 
sented in Section III we use the conventional three-quark 


A. Positive-parity Spectrum 

The results for the nucleon spectrum in the positive- 
parity sector are shown in Fig. [4j Solid horizontal lines 
are added to guide the eye, with their values set by the 
states in basis number 4, since this basis contains all the 
operators studied and has the largest span. 

Of particular interest is the robustness of the varia¬ 
tional techniques employed. While changing bases may 
effect whether or not a particular state is seen, the energy 
of the extracted states is consistent across the different 
bases, even though they contain qualitatively different 
operators. 

Despite the use of 5-quark operators, no state near the 
non-interacting P-wave Ntt scattering threshold is ob¬ 
served. This is understood by noting that none of our 
operators have a source of the back-to-back relative mo¬ 
mentum between the nucleon and pion necessary to ob¬ 
serve an energy level in the region of this scattering state. 

The corresponding eigenvector components for the 
positive-parity states are shown in Fig. [5] as a function 
of basis and variational parameter dt , with to = 17 fixed. 
The values of dt range from 1 through 4. The upper limit 
of dt = 4 was chosen as the largest value for which the 
variational analysis converged for each of the seven bases. 
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n s = 35 + 200 


P-wave N + 7r 


1 -> X 1 + X 2 

2 -> Xl + X2 + XS 

3 -> Xl + X 2 + Xs 

4 -> Xl + X2 + X5 + Xb 

5 —> Xl + X5 + Xs 

6 -> X 2 + XS + Xs 

7 -> X5 + X5 
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FIG. 4. (Colour online). The positive-parity nucleon spec¬ 
trum with various operator bases constructed with 35 and 
200 sweeps of smearing. Horizontal solid lines are present to 
guide the eye and are drawn from the central value of the 
states in basis 4, while the dashed line marks the position of 
the non-interacting P-wave Nn scattering threshold 


The ground-state nucleon is observed in every basis 
regardless of the absence or presence of a particular op¬ 
erator. If xi is present then this provides the dominant 
contribution, with Xs coupling strongly to the ground 
state in bases where Xl is absent. An interesting inter¬ 
play between 35 and 200 sweep smeared xi is observed 
with the smaller source diminishing in importance as dt 
is increased. This may be associated with the Euclidean 
time evolution of highly excited states which are sup¬ 
pressed with increasing dt. 

Turning our attention to state 2, we see that Xi plays 
a critical role in the extraction of the first excited state, 
which is associated with a radial excitation of the ground 
state m- Here the 35 and 200 sweep Xi interpolators 
enter with similar strength but opposite signs, setting up 
the node structure of a radial excitation, xi dominates 
the construction of the optimised operator for this state 
for bases 1 through 5, whereas basis 6 and 7 which lack 
Xi do not observe this state. 

The eigenvectors for state 3, the second excited state, 
are dominated by x '2 components with the same sign 
when this operator is present (bases 1-4,6). This state 
is not observed in basis 5 (where X 2 is absent). Inter¬ 
estingly, in basis 7 which only contains five-quark opera¬ 
tors it appears that it is possible to form this state using 
X 5 components at two different smearings with opposite 
sign. 

We observe that the overall structure of the eigenvec¬ 
tors for each of the three states is highly consistent across 
different bases and different values of the variational pa¬ 
rameter dt. The structure of the eigenvectors can be con- 
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Basis Number/Variational Parameters 


State 2 



Basis Number/Variational Parameters 


State 3 



Basis Number/Variational Parameters 


FIG. 5. (Colour online). Eigenvector components correspond¬ 
ing to the low-lying positive-parity nucleon states. State 1 
corresponds to the ground state, with states 2 and 3 corre¬ 
sponding to the first and second excited states respectively. 
The column numbers denote basis number while the minor x 
axis ticks correspond to the values of the variational parame¬ 
ter dt which runs from 1 through to 4. to = 17 has been used 
throughout. The subscripts 35 and 200 in the legend refer to 
the number of smearing sweeps applied. 
















































































sidered to be a signature or fingerprint of the extracted 
state, and this consistency across bases confirms that it 
is the same state being identified. 

It is fascinating to see that for state 1 in bases 6 and 
7, where X 5 takes the role of the absent Xi operator, 
the values of the two dominant eigenvector components 
(which indicate the mixture of the two different smearing 
levels used) are extremely similar to the xi components 
in bases 1-5. Interestingly, at dt = 2 the error bars for 
the dominant components of states 2 and 3 blow up. As 
we shall explain below, this is due to an accidental de¬ 
generacy in the eigenmasses for this choice of variational 
parameters. 

In order to further test the robustness of our varia¬ 
tional method we conduct a comparison of the masses 
obtained from fitting the eigenstate-projected correlators 
as a function of the variational parameters for each ba¬ 
sis. These results are presented in Fig. [ 6 ] Also shown for 
comparison are the eigenmasses, m a , that result from 
solving the generalised eigenvalue equation of Eqs. |7| 
or ([s]) with c Q = 

Studying state 1, the nucleon ground state, we observe 
that the masses obtained from projected correlator fits 
are approximately invariant across different bases and 
choices of the variational parameter. In contrast, the 
eigenmass lies well above the fitted mass, dropping in 
value as the variational parameter dt is varied from 1 to 
4. While the eigenmass is directly related to the princi¬ 
pal correlator and thus should approach the ground state 
mass in the large time limit, it is clear that the values 
of dt we examine here are insufficient for this to occur. 
It is worth noting that, in bases 6 and 7 where Xi is ab¬ 
sent we see that the eigenmass value rises significantly. 
Nevertheless, the fitted mass remains remarkably consis¬ 
tent with the values obtained in bases 1-5. We emphasize 
how strong the variational parameter dependence of the 
eigenmass contrasts the more consistent structure of the 
eigenvectors. Insensitivity of the eigenvectors to the vari¬ 
ational parameters is a key component of the invariance 
of the masses obtained from the projected correlator. 

Turning to state 2, we see that the eigenmass shows 
similar behaviour to state 1 , lying above the extracted 
mass and dropping with dt. Interestingly, for state 3 in 
bases 1-4 and 6 the eigenmass shows constant behaviour 
for dt = 2 — 4 but systematically lies below the extracted 
mass. In basis 7, the state 3 eigenmass is very different to 
the previous bases, lying above the extracted mass and 
showing a similar downward trend to states 1 and 2 as 
dt varies. 

As for state 1, the fitted masses for states 2 and 3 
provide highly consistent values and uncertainties across 
the different bases and values of dt , with the notable ex¬ 
ception of dt = 2. As observed previously in Fig. [5j we 
see in Fig. [ 6 ] considerably larger error bars at the varia¬ 
tional parameter set {to,dt) = (17,2) in both the eigen¬ 
vector components and projected mass fits for the first 
and second excited states. To understand this, we turn 
to Fig. [TJ where the eigenmasses for states 2 and 3 are 
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FIG. 6 . (Colour online). Comparisons of eigenmasses to 
masses obtained from a projected correlator fit for low-lying 
states in the positive-parity nucleon channel. The column 
numbers denote basis number while the minor x axis ticks cor¬ 
respond to the values of the variational parameter dt = 1... 4. 
to = 17 has been used throughout. The line denoting the ex¬ 
tracted mass is set using basis 4 with dt = 3. 
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FIG. 7. (Colour online). A plot showing the eigenmasses for 
both states 2 and 3, illustrating the accidental degeneracy at 
dt = 2. 


plotted against the variational parameter dt in each basis. 

Note that at (to,dt) = (17,2) there is an approximate 
degeneracy in the eigenmass for states 2 and 3. As a con¬ 
sequence, the corresponding eigenvectors can therefore be 
arbitrarily rotated within the state 2/state 3 subspace 
while remaining a solution to the eigenvalue problem. 
When constructing the jackknife sub-ensembles to calcu¬ 
late the error in the fitted energy, we need to solve for the 
eigenvectors on each sub-ensemble. Due to the approx¬ 
imate degeneracy, the particular linear combination of 
state 2 and state 3 that we obtain for each sub-ensemble 
can vary. Indeed, we observe that the dot-product be¬ 
tween the ensemble average and sub-ensemble can drop 
significantly for dt = 2 in comparison to other values 
of dt. This causes a large variation in the sub-ensemble 
eigenvector components and a correspondingly large er¬ 
ror bar. The simplest way to avoid the problem of this 
accidental degeneracy is to select a different value of the 
variational parameter. 


B. Negative-parity Spectrum 


The negative-parity nucleon spectrum is presented in 
Fig. § Solid horizontal lines have been added to guide 
the eye, with their values set by the states in the largest 
basis (number 4). Once again, while changing bases ef¬ 
fects whether or not we observe a given state, the ex¬ 
tracted states display an impressive level of consistency 
across the different bases. 

The dashed line indicating the energy of the non¬ 
interacting (infinite-volume) scattering-state threshold is 
also indicated with the caution that mixing with nearby 
states in the finite volume can alter the threshold posi¬ 
tion [SUES]- We note here that all scattering thresh¬ 


olds discussed in this section and the next, refer to the 
non-interacting threshold. In contrast to the positive- 
parity results, we do observe a state near the S-wave Nn 
scattering threshold in the negative-parity channel (bases 
5,6,7), also noting that the P-wave Nirn thresholds lie in 
the region of state 3 seen in bases 3, 4 and 5. It is impor¬ 
tant to note that even after the introduction of operators 
that permit access to a state near the low-lying scattering 
state, the energies of the higher states in the spectrum 
are consistent, demonstrating the robustness of the vari¬ 
ational techniques employed. 


n s = 35 + 200 


— S-wave N + tt 

1 -> XI + X2 

2 -> Xi + Xi + X5 

3 -> XI + Xi + Xs 

4 -*• Xi + Xi + Xs + xi 

5 -*• Xi + Xs + Xs 

6 -> X 2 + XS + Xs 

7 -> Xs + Xs 
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FIG. 8. (Colour online). The negative-parity nucleon spec¬ 
trum with various operator bases using 35 and 200 sweeps of 
smearing. Solid horizontal lines are present to guide the eye 
and are drawn from the central value of the states in basis 4, 
since this basis is the largest. The dashed line marks the po¬ 
sition of the non-interacting S-wave Nn scattering threshold. 
The variational parameters used herein are (to,dt) = (17,3). 

Plots of the corresponding eigenvectors for the low- 
lying negative-parity states as a function of basis and 
variational parameter dt = 1... 3 are presented in Fig¬ 
ures [9] and |TU} The upper limit of dt = 3 was chosen as 
the largest value for which the variational analysis con¬ 
verged for all seven bases. The eigenvector components 
for state 0 (when it is present) are dominated by the 
multi-particle operators %5 and Xs> suggesting that this 
state should be identified as a scattering state. The ex¬ 
tracted energy for this state is in the region of the non¬ 
interacting S-wave Nn scattering threshold (which lies 
below the first negative-parity resonant state). The un¬ 
certainty in bases 6 and 7 are relatively large compared 
to basis 5, indicating that the presence of xi may also be 
required to cleanly isolate this scattering state. Indeed, 
we note that in basis 5 there is a significant contribution 
to state 0 from the Xi(n s = 200) operator. 

It is also important to note that either X 5 or x '5 can 
be the dominant interpolator exciting this lowest-lying 
state. Given that X 2 is predominantly associated with the 
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Basis Number/Variational Parameters 

FIG. 9. (Colour online). Eigenvector components correspond¬ 
ing to State 0 which is in the region of the non-interacting 
S-wave N + 7r scattering threshold. The column numbers de¬ 
note basis number while the minor x axis ticks correspond to 
the values of the variational parameter dt which runs from 
1 through to 3. to = 17 has been used throughout. The 
subscripts 35 and 200 in the legend refer to the number of 
smearing sweeps applied. 


third state in the positive-parity sector at 2.4 GeV one 
might naively expect X 5 would be associated with S-wave 
scattering states near 2.7 GeV. Remarkably it creates 
a scattering state near 1.35 GeV. Thus one should use 
caution in predicting the spectral overlap of five-quark 
operators by examining the spectral overlap of the pion 
and nucleon components of the five-quark operators sep¬ 
arately. In light of the quark field operator contractions 
required in calculating the full two-point function this 
result is not surprising. 

In accord with previous studies nra, we find that 
the Xi interpolating field is crucial for extracting state 
1 , associated with the lowest-lying negative-parity reso¬ 
nance, as we do not observe this state when xi is absent 
as in bases 6 and 7. As expected, Xi provides the dom¬ 
inant contribution to state 1 , which is associated with 
the 5ii(1535) in Nature. Similarly, we see that %2 has a 
high overlap with state 2, the next resonant state. Basis 
5 does not see state 2 due to the absence of X 2 - However, 
unlike state 1 , there is an important mixing of xi and x '2 
in isolating the eigenstate. It is interesting to note that 
in basis 7 we are able to form this state by combining xs 
and X 5 - 

The consistency of the eigenvector structure for the 
low-lying states 1 and 2 is strong. Despite the appear¬ 
ance of a state near the S-wave Ntt threshold, state 0 in 
basis 5, the eigenvector components for state 1 are re¬ 
markably consistent with those in other bases where this 
lower-lying state is absent. If we look at basis 6 , where 
state 0 is present but state 1 is absent, the eigenvec¬ 
tor components for state 2 are in good agreement with 



Basis Number/Variational Parameters 


State 2 



Basis Number/Variational Parameters 


State 3 



Basis Number/Variational Parameters 

FIG. 10. (Colour online). Eigenvector components corre¬ 
sponding to low-lying negative-parity nucleon states. States 
1 and 2 correspond to the two lowest-lying resonant states, 
while state 3 interestingly lies in the region of the P-wave scat¬ 
tering thresholds. The column numbers denote basis num¬ 
ber while the minor x axis ticks correspond to the values of 
the variational parameter dt which runs from 1 through to 3. 
to = 17 has been used throughout. The subscripts 35 and 200 
in the legend refer to the number of smearing sweeps applied. 
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those from other bases where the lower-lying state 0 is 
not observed. This demonstrates that, with a judiciously 
chosen variational technique, a reliable analysis of higher 
states in the spectrum can be performed even if states 
associated with the low-lying scattering states are not 
extracted by the correlation matrix analysis. 
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FIG. 11. (Colour online). Comparisons of eigenmasses to 
masses obtained from a projected correlator fit for state 0, 
which is in the region of the non-interacting S-wave Nn scat¬ 
tering threshold. The column numbers denote basis number 
while the minor x axis ticks correspond to the values of the 
variational parameter dt = 1... 3. to = 17 has been used 
throughout. The line denoting the extracted mass used has 
been set using basis 5 with dt = 3. 

State 3, which lies in the region of the non-interacting 
P-wave Nttit scattering states in the channel, also shows 
good agreement across bases and variational parameters. 
The dominant eigenvector components show that this 
state is formed from a mix of Xg and Xi- If is worth not¬ 
ing that very early choices of the variational parameters 
such as (to, dt) = (17,1) do not allow sufficient Euclidean 
time evolution to cleanly isolate this state. The correla¬ 
tion matrix has more states participating in the analysis 
than the dimension of the basis leading to contamination 
from unwanted states and hence spurious results. The 
different structure for the state 3 eigenvectors at these 
early variational parameter sets illustrates the need to 
allow sufficient Euclidean time evolution to occur. 

The comparison of the fitted masses as a function of 
variational parameter dt across the different bases for 
the negative-parity sector is shown in Figures [TT] and [l2| 
Again, the eigenmasses are plotted for comparison. As 
before, we observe for all the states the fitted masses are 
consistent across the different bases and values of dt. In 
contrast, the eigenmasses for the negative-parity states 
all show some variation with dt to different extents, with 
the values typically lying well above the extracted ener¬ 
gies. 

Finally, we observe that whenever X5 is present, ei- 
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FIG. 12. (Colour online). Comparisons of eigenmasses to 
masses obtained from a projected correlator fit for low lying 
states in the negative-parity nucleon channel. The column 
numbers denote basis number while the minor x axis ticks 
correspond to the values of the variational parameter dt = 
1... 3. to = 17 has been used throughout. The line denoting 
the extracted mass used has been set using basis 4 with dt = 3. 
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ther a state near the S-wave Nn scattering threshold, or 
a state lying in the region of the P-wave Nnn scatter¬ 
ing thresholds is extracted. This indicates the presence 
of the vector di-quark in the interpolator may play an 
important role in scattering state excitation. It is per¬ 
haps surprising that basis 4 fails to see a state near the 
lowest-lying scattering threshold in the sector, despite 
being the largest basis. We believe this is due to the 
spectral strength available to the scattering state being 
relatively low. The overlap of the scattering state with 
the operators is not high enough to compete with the 
large spectral strength imparted to the low-lying resonant 
states when both \i and X 2 are present. We note that 
the only time our local (three-quark or five-quark) op¬ 
erators overlap with a meson-baryon state is when both 
hadrons are at the origin. The probability of this occur¬ 
ring is proportional to 1/V 1 2 3 4 5 6 . After taking into account 
the spatial sum in Eq. [lj this results in a \/V suppres¬ 
sion of multi-particle states in the correlator amplitude 
Q(t) [ 0 ]. Indeed, it seems to be relatively difficult to 
extract a state near the S-wave Nn state with our lo¬ 
cal five-quark operators, suggesting that scattering state 
excitation is best achieved by explicitly projecting the 
momentum of interest onto each hadron present in the 
scattering state. 

VII. CONCLUSIONS 

We have investigated the role of local multi-particle in¬ 
terpolators in calculating the nucleon spectrum by exam¬ 
ining a variety of different bases both with and without 
five-quark operators. 

The variational techniques herein employed, demon¬ 
strate that fitting a single-state ansatz to optimised 
eigenstate-projected correlators provides a method to re¬ 
liably extract energies in both the positive and negative- 
parity channels. While the selection of states that are 
observed varied between bases, when a given state is seen 
the extracted energy agrees across qualitatively different 
bases. 

Furthermore, the structure of the eigenvector compo¬ 
nents and the corresponding fitted energies for the states 
observed are shown to be highly consistent across dif¬ 


ferent bases and choices of the variational parameters, 
despite the markedly different interpolators used in the 
various bases. We found that an approximate accidental 
degeneracy in the eigenmass at (to, dt) = (17,2) for states 
2 and 3 led to a large increase in the uncertainties for the 
corresponding energies and eigenvector components. 

While we did not observe any positive-parity scatter¬ 
ing states, in the negative-parity sector we found that X 5 
was crucial to obtaining an energy in the region of the 
non-interacting S-wave Nn. Even with the use of local 
five-quark interpolators the uncertainties on this thresh¬ 
old state were relatively large compared to those of higher 
states. Future studies will include multi-particle opera¬ 
tors with explicitly projected single-hadron momenta in 
the variational basis to facilitate better excitation of scat¬ 
tering states, including those in the positive-parity sec¬ 
tor. An interesting feature of our negative-parity results 
is that the energies of the extracted states are consistent 
across all bases in which the state is observed, regardless 
of the presence (or not) of a state in the region of the 
lower-lying non-interacting scattering threshold. This 
suggests that by using the techniques described herein, 
one does not need to have access to the aforementioned 
low-lying states to reliably extract energies closely related 
to the resonances of Nature. 
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